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NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

TECHNICAL MEMORANDUM X-JOO 

AN ANALYTICAL INVESTIGATION OF ABLATION* 

By Bernard Rasnis and Russell N. Hopko 


SUMMARY 


An analytical procedure is described which enables the derivation of 
effective heat of ablation relationships for any type of boundary layer 
from transpiration cooling results. The procedure enables the inclusion 
of such effects as the ratio of wall temperature to local stream tempera- 
ture, surface radiation, and surface combustion. 

The predicted effective heats of ablation for a three-dimensional 
laminar stagnation boundary layer for Teflon material were in agreement 
with those derived from tests conducted at boundary-layer enthalpy 
potentials of 800 and approximately 7,000 Btu/lb. 

The predicted equilibrium surface temperatures on nonablating surfaces 
behind an ablating material were in agreement with the values derived 
from tests conducted with Inconel cylinders having Teflon hemispherical 
nose pieces . 


INTRODUCTION 


Previous experimental investigations (refs. 1, 2, and 3) have indi- 
cated that ablating materials have considerable promise as a cooling or 
heat-blocking system for high-speed vehicles, from a standpoint of sim- 
plicity as well as efficiency. Reference 4 gives the framework of an 
excellent theoretical treatment which shows how the similarity parameters 
of ablation are related to the parameters of transpiration cooling. The 
analysis of reference 4 does not, however, include the effects of dif- 
ferent values of the ratio of wall temperature to local stream temperature 
nor does it include ablating materials having specific heat values dif- 
ferent from those of air. 

The purpose of this paper is to show how expressions relating the 
ablation parameters to the material and boundary- layer characteristics 
may be derived for any type of boundary layer, simply and directly from 
transpiration cooling results. 

* Title, Unclassified. 
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SYMBOLS 


C 

C P 

F 

f 

G 


H 

h 

heff 


L 

M 

m 

%u 

N pr 

%e 

N St 

q. 


velocity gradient 
specific heat, Btu/lb-°F 


ratio of coolant mass -flow rate to local air mass -flow rate 
fraction of vaporization 

average value of coolant flow rate, lb/(sq ft) (sec) over 
injection area 

enthalpy, Btu/lb 

heat-transfer coefficient, Btu/(sq ft)(sec)(°F) 
effective heat of ablation, Btu/lb 


vap 


'^ c p,melt (^vap 



latent heat, Btu/lb 
molecular weight 

mass ejection rate, ib/(sq ft) (sec) 
Nusselt number 


Prandtl number 
Reynolds number 
Stanton number 

heat flux, Btu/(sq ft) (sec) 
heat transfer from gas layer to 


ablating surface (see eq. (l) 
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S surface area, sq ft 

T temperature , °F 

t wall thickness, ft 

V velocity, ft/sec 

W mass fraction of foreign gas 

w mass flow, lb/sec 

x porous length plus downstream length, ft 

x Q porous length, ft 


A 

e 

T 

^comb 


heat of combustion, Btu/lb 
emissivity 

boundary -layer recovery factor 

^Icomb 

^aero 


^rad 

^rad ~~ q 

uaero 

v kinematic viscosity, ft^/sec 

p density, lb/cu ft 

a Boltzman's constant 

t time 

Subscripts : 

aero aerodynamic 

air air 

aw adiabatic wall 

c coolant 
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comb combustion 

e equilibrium 

gas gas 

Z local 

Ze portion of leading-edge area that is supplying coolant 

m mixture 

melt liquid film or melting 

o zero coolant flow or no ablation 

p porous area 

rad radiation 

solid solid material 

sub sublimation 

t total 

vap vaporized layer or vaporization 

w wall 

w,air air at wall temperature 

£ downstream of porous area 

co free stream 

ANALYSIS 

In figure l(a), there is shown a schematic diagram of a body under- 
going ablation at the stagnation point. The body is heated and melts, the 
melt forming a liquid film. This film or a portion thereof is then vapor- 
ized, the vaporized material forming a gas layer. This process is herein- 
after referred to as ablation by melting and vaporization. In figure l(b), 
typical temperature variations are shown for the different phases . 
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Defining q^ as the heat transfer from the gas layer to the ablating surface, the general 
-balance equation for the body of figure 1 can be written (ref. 5) as 


^ " “Vap L vap + 3^, melt (tap “ ^melt) + “melt j c p,melt ^-vap " T melt) + ^melt + °p, solid (^melt " T soli4) W 

\ r here ^elt i s °i material melted per second and iyap is the amount vaporized per 


second. 


If the heat transfer to a nonablating surface which is at the same temperature as the 
iblating surface is denoted by <la er0 j equation (l) can be written as 


|“vap Wap + ^p^melt ( T vap ~ T meltJ] + “melt ^p^melt^ 


^vap “ T melt) + ^elt + c p 


, solid ('^mi 


.elt " T solid j 


For a melting-vaporizing system if f is defined as that fraction of the melt which is 
vaporized, or 


^vap ~ -^^melt 


equation (2) becomes 


-JL 


“melt 


^ |^vap + ^p^melt (^vap “ ^meit)j + ^p^melt (^vap " ^melt) + ^melt + solid (■ 


t T melt “ T solidJ 


If the term that is multiplied by m^p^ ( i en °i e( i by ^ov an( ^ 311 effective heat absorp- 
tion parameter that takes into account both the heat absorbed by the ablating material and the 
heat shielding effect of its vapor is defined as 


heff = 

%ielt 
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there results 


^ _ ^lelt^Qv _ ^ov 

^ero ^aero ^eff 


Since the term 


is analogous to the shielding-effect ratio 


N 


St 


N, 


St, aero 


^ero 

obtained from transpiration cooling, then 


N, 


St 


A ov 


%t,aero h eff 


(4) 


The transpiration cooling flow rate correlation parameter defined as 


N, 


c P,c 

St, aero °p , Z 


m vap 

(pV)i 


J);C 


%.ero 


°P ,7- 


(p^)l c p,z(^’aw “ 


can be written as 


F Cp,c “vap „ ^ c P,c 

ivr “ “ V^aw " ^ ,&±rl ~ 

%t,aero °ji,l < 3aero' 


aw_ ~ ^w.air) c P,c 
^ero °p , l 


~ ^w,air) °p,c 
^ff °V) 1 


but from equation (4) 
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N, 


h eff = h ov 


St, aero 


N, 


St 


Thus, 

F °p,c _ ffoaw ~ ^w,air) %t °p,c 
*%t,aero °~p,l ^ov %t,aero ^p,! 


and 


F °P,c 
Nst,aero c p,Z 

N St 

%t,aero 


^aw ~ ^w,air) c p,c 
hov , I 


( 5 ) 


Hence, once the relationship between Ngt/%t,aero an< 3- 
^F/Ng-(-^ aero ^ ^, 0 /^, 1 ) known from transpiration data, the transfor- 
mation relationships (eqs. (4) and ( 5 )) enable the solution of h eff in 
terms of the enthalpy difference H aw - H v a ^ r and the intrinsic heat 
capacity of the material h ov . Thus, once the shielding relationship is 

known for any boundary -layer type and vapor, evaluation of the material 
overall effectiveness depends only on the material properties in the 
solid and molten states, and the quantity f . 


A discussion of f (percent of vaporization) is given in reference 6 
Briefly, f increases with flight speed at constant pressure and 
decreases with increase in pressure at constant flight speed. For any 
given flight conditions, f increases with the Prandtl number of the melt 


For the case where radiation, combustion, and other effects modify 
the reference heat-transfer rate, an actual effective heat of ablation 
parameter is defined as 


h eff 


* 


h eff 

1 + ill + t ) 2 + . • ■ % 


where 


= + . 


Aq 


Qaero 
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the plus or minus sign depending on whether Aq adds or subtracts heat 
from ‘laero * 

For example, the effect of radiation from the surface is to 
decrease the reference heat-transfer rate; thus, 


Aq ra d 

^rad ~ 1 " „ 1 

^aero 


0.48 

tW T * 'll 

1 4 


\1000jj 



%.ero 


and is equal to 1 when Aq ra ^ is zero (e 1 /^ = o) and is equal to 
zero when the numerical value of the parameter e^/^T v is equal to the 

radiation equilibrium temperature. Some representative values of 
1 - q ra £ are shown in figure 2. 


The actual mass losses m 


.* 


m = 


are thus defined as 
%,ero 




and the ratio 


m 

m 


^ff 

heff" 


The effect of combustion would be to increase the reference heat- 


transfer rate; thus, 


1 + 


^comb 


= 1 + 


^^comb 


^aero 

The following expression has been derived in reference 7 : 


al. 


^aw ” ^v,air + ^ 

®aw _ ^w^air 


( 6 ) 
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heat transfer from gas layer to ablating surface when combustion 
occurs 

ratio of molecular weight of product to that of oxidizer 

mass fraction of oxidizer evaluated at outer edge of boundary 
layer 

heat of combustion per pound of product 

The above expression is valid only for laminar boundary layers and 
assumes a Lewis number of 1.0 (Prandtl number = Schmidt number). It 
applies to the cases where the mass flow of oxygen to the surface is a 
maximum. This maximum flow rate of oxygen is assumed to be completely 
consumed in the reaction. Below some small value of the mass-flow rate 
of the combustible coolant corresponding to stoichiometric surface com- 
bustion, the expression is no longer valid. 

Upon rearrangement , equation (6) can be written as 


where 

K 


A 


V ~ ^ = Aqcomb 
%.ero ^.ero 


kpy 

^ff 


KW 


for, air 





Since 

_ fr-QV 

^ero ^eff 


and 


thus, 


^w,air 



ov 


KW 


^ + ^comb ” ^ + 


A ov 


y^airy 


^eff ( ~^aw " ^,air\ °p,c 
^ov / °p , l 
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and is equal to 1.0 for no combustion (A = 0) and approaches maximum 
values for stoichiometric combustion at the surface. Some representative 

KA 

values of 1 + Tj com b are shown in figure 3. The curve for - 15 

*V,air 


represents the chemical reaction 


2H 2 + 



2H 2 0 


The curve for 


KA 


v 


= 20 represents the combustion of Teflon. Refer-* 


air 


ence 6 gives the reactions and A values for this material. 


Derivation of Shielding -Effect Relationships 

In the preceding section, it was shown that effective heat of 
ablation relationships may be obtained directly from transpiration cooling 
results. In this section, several available theoretical and experimental 
transpiration cooling results for several types of boundary -layer flows 
will be transformed into the pertinent ablation cooling parameters by 
means of equations (4) and (5). 

Two-dimensional turbulent .- Reference 8 gives experimental results 
for helium-to-air and nitrogen-to-air transpiration cooling of a two- 
dimensional 20° half -angle wedge for a free -stream Mach number of 2.0 and 
stagnation temperatures ranging from 1,295° F to 2,910° F. The results 
of reference 8 are presented in that paper as the shielding-effect ratio 

N of jr 

— — as a function of the flow-rate parameter . However, 

^St,aero ^St,aero 

it can be shown that to a good approximation (see appendix) the coolant 

weight requirements, for the same cooling ratio, or 


N- 


St 




NSt_\ 


St,aero/li e iium_t 0 - a ir ' ^' aero ^air-to-air 


are inversely proportional to the specific heats of the coolants, or 



helium 


K,= 


)air 
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Thus, the data of reference 8 were modified by multiplying the 


parameter 


F 

%t,aero 


by the corresponding ratios of 


C PjC 
C P ,1 


where 


c V,c } 


the coolant specific heat, was evaluated at the cooled wall temperature, 
and c^ 7 , the stream specific heat, was evaluated at the local tempera - 

ture . The modified results are shown in figure 4. 


In reference 9> there is derived a simple expression for the 
shielding effect of transpiration for a two-dimensional turbulent boundary 
layer. The expression can be written as 


-p,c 


N, 


St 


N. 


N, 


St, aero 



St, aero °p ,2 
F c Pj c > 
St, aero c p,’i / 


The curve computed from this expression is also shown in figure 4 
(dashed curve). The theoretical values of reference 9 are somewhat 
higher than the experimental data of reference 8 . It should be noted, 
however, that the theory of reference 9 is for T w /T ^ equal to 1.0, 

whereas the data of reference 8 were for values of T w /T£ ranging from 

approximately l /8 to slightly less than 1 . 0 . Since the shielding -effect 
ratio would decrease somewhat for decreasing T v /T^, the difference 

between the two curves of figure 4 is probably due to this effect. 
Unfortunately, the experimental scatter in the data of reference 8 pre- 
cludes any delineation of the T w /T ^ effect. 


Two-dimensional laminar .- Reference 10 gives theoretical transpira- 
tion cooling results for air-to-air and values of T w /T^ of 1 . 0 , l/ 2 , 

and l/4. The results of reference 10 are presented in terms of a heat- 


N- 


transfer parameter 


Nu,w 


/ 


and a nondime ns ional flow rate outward from 


N- 


Re,w 


the porous surface f w . Although reference 10 computes the cases for 
values of 


r w/ T ! 


other than 1.0 by the procedure described it was 


found that by dividing 
the ratio 


L w 


by the ratio 


( T w/ T 0 


0.095 


the values of 
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( N Nu,v/y N Re,v) fw ^ Q 

( N Nu,w// N Re,w ) fw=0 


all fell on a single curve as shown in figure 5 . 


Since the preceding ratio is analogous to the ratio 


»st 

%t,aero 


and 


v N St,aero / t/T z =1.0 


(~^v)Npr,w 


N- 


Nu,v 


Ni 


R e,w/ f =0 
w j 


T w /Tj=1.0 


the results of reference 9 vhen plotted as the ratio 
the ratio 


N St 

aero 


against 


r , F A C P,e 

^St, aero lt0 °V,l 

( V T i) 0,095 


provides the shielding curve for a two-dimensional laminar boundary layer 
applicable also to the cases where is different than 1 and also 

where the transpired coolant has properties different than air, . 

Two- and three-dimensional laminar stagnation ,- Reference 11 gives 
theoretical transpiration cooling results for air-to-air and V T * 

equal to 1.0, 0.5, and O. 25 . The results of reference 11 are given in 
terms of the heat-transfer parameter N^ u w //^Re w ancl a non ^^ ens i° na l 

normal velocity outward from the porous surface V w ^y^C. However, for 
the two-dimensional stagnation case 


( _f w) 
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and. for the three-dimensional stagnation case 


(”^w) 



As shown in figure 6 , the results of reference 10 when plotted as 
the ratio 


N Nu, v / / (%u,v/^) v 


against the ratio 


Vw/AvO 


( V T 2 )' 


0.19 


reproduce as single curves. Since 


N, 


St , aero/ ^ /T ^=1 . 0 


•w 


Nr 


vTW ,v 


( n Nu,w ^ 
\\/ N Re,wy 


V T 2 =1 -° 


the results of reference 10 when the ratio 


function of the ratio 


N, 


St 


N, 


St , aero 


is plotted as a 


^c 


\ N St,aero/m / T 


provide the shielding curves for the two- and three-dimensional laminar 
stagnation boundary -layer cases, applicable to all coolants and values 
of T w /Tj. Use of the transformation equations (4) and ( 5 ) converts the 

shielding curves into the curves applicable to ablating materials. 
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Figure 7 shows the effective heat of ablation curves that were 
obtained from the shielding curves previously discussed for the case 
where T w /T 1 is 1. 

Comparison of Experiment and Theory 
Figure 8 shows h , e ff of Teflon as a function of the ratio 


(^aw - Hw ,air) °p,c 

(t w /t z ) 0 - 19 Cp,Z 

for three-dimensional laminar stagnation flow. The circular symbols 
denote the experimental values derived in reference 3. The enthalpy 
potential ratio was evaluated using a value of 900 Btu/lb for h ov . 

The values of T v and CpjC ve re obtained from figures 9 and 10, respec- 
tively. The curves of figures 9 and 10 were obtained from references 12 
and 13, respectively. Emissivity values for Teflon were obtained from 
reference 14. The values of c p ; 2 vere obtained from reference 15. The 

values of T^T^ and Z vere 0.498 and 0.93 for ‘the ethylene - 

jet test and equal to O.I67 and 0,5 for "the arc- jet test. Since the 
products formed by the depolymerization of Teflon are a gas, f is 
equal to 1. In some static tests, a liquid phase has been observed } 
however, in computing the values shown in figure 8, f equal to 1 was 
used. 


The vertical lines shown in figure 8 represent the data of refer- 
ence 16. No values of T w /T ^ and Cpjc/^p,! were available for the 

data of reference l6. However, reference l6 uses (“air/Mgas)^ 4 to 
correlate the results for different material vapors j thus, it was assumed 
that the ratio Cp^ c / C p, Z f(^w/^l) ^ was O.76 ^the value of 

(^air/^o-as) f° r Teflon). Thus, the enthalpy potential values of 
reference l6, H aw - H^ air , were multiplied by O.76 to conform to the 

parameters of figure 8. It should be noted that the data of reference 16, 
as shown in figure 8, would give an intercept value 900 Btu/lb (calculated 
value of h ov for Teflon) instead of 75O Btu/lb as shown in reference l6. 


Downstream Cooling Effect of Ablating Materials 

A problem related to ablation cooling for which solutions are 
obtained by essentially the same procedure as described in the previous 
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sections is that of the downstream cooling effect on a solid surface 
behind an ablating material. In downstream cooling, the coolant is 
ejected into the boundary layer forward of the surface which is to be 
cooled and is swept back over the rearward section by the air flow. In 
the following sections it will be shown how the temperature distribution 
results on a nonporous plate behind a region with transpiration cooling 
can be correlated in terms of the parameters pertaining to ablation 
cooling. 

Two-dimensional laminar compressible flow .- In figure 11 , there are 
shown the results of the analysis of reference 17 for the temperature 
distribution on an insulated nonporous plate behind a region with trans- 
piration cooling. The results are given in terms of the wall cooling 
efficiency parameter (t w - T c ) j ( T a v - T c ) as a function of the distance 
behind a porous region for a number of fixed flow rates. A family of 
curves is obtained by this procedure, each curve pertaining only to a 
certain flow rate. Since x represents the sum of the porous and down- 
stream areas and x 0 represents the porous area, the quantity (x/x 0 - l) 

is equal to Sg/Sp, which is the ratio of the downstream cooled area to 
that of the porous area. 

Since , by definition 


- f O = ^ 


where G is the nondimens ional flow rate averaged over the area of 
injection, m/(pV)j, the term 

G 


x 

*o 


- f Q .. S p/%e,°° G 

- " S l/ S P ~ s i/%^ 


= Ft ,/N 


!y H Re,c 


Dividing by %t,aero,^^Re « (ref. 18) gives the downstream flow-rate 
correlation parameter 


. N< 


St,aeroy 


The single correlation curve obtained by this procedure is shown in fig- 
ure 12 . The values of aero)^ w ^re multiplied by the ratio 

Cp c jc p 1 . This permits application of the results to any coolant. 

(See appendix.) 
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Two-dimensional turbulent compressible flow .- Also shown in fig- 
ure 12 are the faired experimental results of reference 19 . These 
results were derived from tests conducted with a two-dimensional wing 
segment having constant chord thickness. The wing segment was con- 
structed with a semicircular porous leading edge through which helium was 
ejected. The tests were conducted in a Mach number 2.0 ethylene -heated 
high-temperature air jet for stagnation temperatures ranging from 500 ° F 
to 2, LOO 0 F and for angles of attack of 0°, +5°, and tl5°* 

Correlation in terms of ablation cooling . - Since 


°p,c _ °p,c _ G/Sp Sp %t,aero,p c p,c 

St, aero i. ^St,aero,| %t,aero,p %t,aero,| C P^ 


then 


f F \ °P,c _ / F \ °p,c S p N St,aero,p 
\%t,aero/| c p,2 \%t,aero/p C p,t %t,aero,£ 

where the factor ^F/Ngt >aero ) is the flow-rate correlation parameter 
for the area of injection. Since, 


F _ ^(^aw ~ Ay, air) 


N, 


St, aero 


x eff 


the downstream flow-rate correlation parameter can be written as 


iN, 


St,aeroy 


f(Haw - Hvf^air) 


heff 


le 


r c p,c|/ s Ze < \ f %t,aero, Ze j 
~ » V vl /\ St, aero, | / 


( 7 ) 


and, thus, directly applies when ablating materials are used at the 
leading edge . 

Since the quantity G is an average value over the injection area, 
the average mass loss rate m is determined from 
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T 3-aero.Ze 

m = ’ — 

n eff,le 

and is an average value for the portion of the leading-edge area which 
is supplying the vaporized mass that cools the downstream surface . The 
coolant temperature T c which is the temperature of the ejected coolant 

at the junction of the injection and downstream areas becomes the vapori- 
zation temperature T va p or the sublimation temperature Tg^ depending 

on whether the material melts and vaporizes or directly sublimes. For 
the case where the material directly sublimes (f = 1.0), the use of fig- 
ure 12 is straightforward. The material bgff at the leading edge is 

calculated for f = 1.0 and air equal to c p :r x T sub* downstream 

flow-rate correlation parameter is then computed for the various downstream 
positions and T| is consequently determined from the curves of figure 12. 
For the case of melting vaporization where all the melt vaporizes at the 
leading edge (f = 1.0), the procedure is the same as for the sublimation 
case, except is now equal to Cp^ T V ap* 

For the case where f is less than 1, some of the melt flows back 
over the afterbody, the length of melt depending on the melt heat capacity 
and the heat input into the melt as it flows along the afterbody. In 
this case, the downstream cooling procedure for f = 1.0 is applied at 
the position where the melt terminates. The afterbody surface that is 
covered by the melt is at a temperature that is essentially equal to the 
temperature of the underside of the melt. 

The circular symbols shown in figure 12 are the experimental results 
obtained from a series of tests conducted in a Mach number ^.0 ceramic- 
heated jet. The test configuration was a 5/8-inch-diameter hemisphere- 

cylinder. The cylindrical section was approximately 2~ inches long and 

constructed from l/l6-inch-thick Inconel. The hemispherical cap was 
Teflon material. The Inconel section was instrumented with thermocouples. 
The details of these tests are given in reference 20. Unfortunately, 
reference 20 presents only the value of m at the stagnation point, and 
does not provide the experimental information to obtain m from equa- 
tion (8) . However, 


fe) 5 - i¥) t 

where w is the total weight loss of the nose material per second, and 
Cp c | , ^o,| > and are the specific heat of the nose material vapor, 
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heat-transfer coefficient, and cooled area of the downstream position. 

The values of w were obtained from Douglas Aircraft Company, Inc.; 
the values of T c ^T c = T w j and Cp,c> from figures 9 and 10; the 

values of T| , from extrapolation of the temperature data of reference 20 
by means of a plot of log T against 1 /t. 

It should be noted that the afterbody surface temperatures shown in 
figure 12 are steady-state or equilibrium values. For the case where the 
afterbody surface is wetted by the coolant but has not as yet reached 
steady-state temperatures, the temperature during the transient period is 
given by the following equation: 



The values of (%t/%t,aero)| are shown in figure 13 for both the 
laminar and turbulent cases as a function of the parameter 


, aercy | 


°P ,c 

c p ,2 


A W - H v ^ a i r \ 
[ heff ) 


7e 


N St,aero , le c p,c 
Ngt,aero,£ c p,Z 


Radiation effects are accounted for by correcting the term Iiq ^ . 

Since 


^ero ” ^rad 


or 


a o (™aw “ T v) “ ^aero(^av ^w) ~ a€ -^w 


there results 


h 


o 


- h 


aero " 


qgT w 4 
(^aw “ T w ) 


CONFIDENTIAL 



CONFIDENTIAL 


19 


or for the downstream surface 


^o^ ^aero,| 



CONCLUSIONS 


The following conclusions are based on the results obtained from an 
analytical investigation of ablation cooling: 

1. Effective heat of ablation relationships for several types of 
boundary layers applicable for all wall-to-local temperature ratios can 
be derived simply and directly from transpiration cooling results, 

2. The predicted effective heats of ablation for a three-dimensional 
laminar stagnation boundary layer for Teflon material are in agreement 
with those derived from tests conducted at boundary-layer enthalpy poten- 
tials of 800 and approximately 7,000 Btu/lb. 

3. The predicted equilibrium surface temperatures on nonablating 
surfaces behind an ablating material are in agreement with the values 
derived from tests conducted with Inconel cylinders having Teflon hemi- 
spherical nose shapes. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va. , April 25, i960. 
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APPENDIX 

DERIVATION OF THE CORRELATION PARAMETER c Pj> c/ c p ; Z 


In the sections dealing with the derivations of the shielding 
expressions , it was indicated that the parameter 


N. 


St, aero 



[where 




the coolant or vapor specific heat, was evaluated at the 


cooled wall temperature and the l° ca l stream specific heat, was 

evaluated at the temperature of the outer edge of the boundary layer^ 
essentially correlated all the cooling results regardless of the type 
of coolant used. 


For the case of foreign gas ejection, the gas layer adjacent to the 
cooled surface is made up of both molecules of foreign gas and the local 
stream. The specific-heat capacity of this mixture will be denoted as 
Cp^ m and is defined 


°p,m “ ^°p,c + " W )°p,w 


where W is the mass fraction of the foreign gas evaluated at the wall 
and c„ w is the local stream specific heat evaluated at the wall 

p, w 

temperature . Thus 

+ (1 - W)^ 

C P ,1 c V,l 

and 

G p> m = 

c p,c c P,c/ c p,l 

Reference 21 gives theoretical results for W as a function of f w 
for hydrogen-to-air cooling for both the two- and three-dimensional 
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laminar stagnation cases. Reference 22 gives results for the hydrogen- 
to-air cooling of a two -dimensional laminar layer. Reference 23 gives 
results for the helium-to-air cooling of a two-dimensional laminar layer 
The results of these references are given in figure 14 where c p,m/ c p,c 

is plotted against Fj Ng^ aero (fw is analogous to this parameter^. 

The results indicate that for the four cases shown, Cp m /Cp^ G becomes 
equal to 1.0 when F^ Ng^ aero g rea ^er than 0.7* Thus, since 



^St, aero ^eff 


the generalization of the ablation formulas to all types of vapors, 
through use of the factor Cp c Jc p ^ , is probably very good when 

-^(^aw " -^w,air) ^ ^eff 


For steady-state or equilibrium conditions (radiation and conduction 
assumed negligible), the heat balance equation for transpiration cooling 
is 


N St = F C P, m ( T w - ?c) 

^St,aero ^St,aero c p,Z (T e - T w ) 

or in view of the preceding argument 

%t F °P,C (T„ - T c ) 

%t,aero ^St,aero c p , Z (^e “ ^w) 


where T e is the boundary -layer recovery temperature for the mixture of 
coolant and stream flows and is defined as 


■^e T c + (^t " ^c) 


The term q e 


is the temperature recovery factor and is defined as 
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reference 9 gives experimentally determined values of Tj e for both air- 

to-air and helium-to-air transpiration cooling of a two-dimensional 
turbulent boundary layer as a function of the flow-rate parameter 
F /%t,aero* The resu lbs reference 9 were modified by multiplying 

F'/Nst aero by the ratio G p^ c / C p , Z anci Plotting T] e as a function of 
(F/Ngt aero) ( c p cj ° p l) * This plot is shown in figure 15 and it is 
seen that rj e is essentially linear with ^F^Ngj. aero ) ^ c p cy/ c p i) an( ^ 
can be expressed by the relationship 


^e 


^aw 


i 


0.015 


N, 


St, aero y 



where i] aw is the temperature recovery factor for no mass transfer. 

Thus, for the same T w and T c (coolant storage temperature), the 
ratio N St /N st ^ aero is only a function of (F/N st ^ aero ) (cp, c /Cp, z ) since 
T e is also specified by this ratio. 

In correlating the experimental results of references 8, 9; and 19, 
G p,Z was evaluated for the temperature at the outer edge of the boundary 

layer. At high temperatures, however, the Cp j varies considerably 

with small changes in temperature and pressure. In order that the cal- 
culations be not incorrectly influenced by these variations, the Cp ^ 

which was used at high values of enthalpy potential was the average 
specific heat for the temperature range through the boundary layer. This 
average value was calculated from 


H 7 


-V,l 


Hw 


T l - T v 
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Figure 2.- Variation of radiation effect with aerodynamic heat flux and surface temperature 
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Figure 5 .- Transpiration shielding results of reference 10 for a two-dimensional laminar 

boundary layer. 
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Figure 6.- Transpiration shielding results of reference 11 for two- and three-dimensional 

stagnation laminar boundary layers. 
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Figure 7*- Effective heats of ablation as a function of the ablating material and 

boundary-layer characteristics Tw/T i = 1. 
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0 Experimental (ref. 5) 

! Experimental (ref. 16) 
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Figure 8.- Comparison of theory with experiment for Teflon. 
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Figure 9>- Variation of surface temperature with ablation rate for 

Teflon. Reference 12. 
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Figure 11.- Temperature distribution on an insulated nonporous plate behind a region with 

transpiration cooling. Reference 17- 


CONFIDENTIAL 



aw 



Figure 12.- Variation of downstream cooling as function of upstream ejection and cooled area 
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Figure 13- - Variation of heat reduction 


Faired data of ref. 19* 
two-dimensional turbulent 
Meo = 1.5 to 2.5 


Ref. IT, two-dimensional laminar 
compressible theory* = 2.0 
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function of upstream ejection and cooled area 
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Figure 14.- Variation of mixture specific heat with flow-rate parameter. 
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Figure 15 *- Variation of boundary -layer recovery factor with the flow-rate correlation parameter. 
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